A. REMARKS 

Claims 1-4, 6-9, 20-23 remain pending. Claims 24-26 are added by this 
response. No new matter is added by the newly presented claims. 

B. Rejections under 35 U.S.C. 112 

Claims 1-4, 6-9 and 20-23 were rejected under 35 U.S.C. 112, second 
paragraph, as being indefinite because the specification does not explain how the axial 
dispersion coefficient would be determined in a given system. 

The Applicants respectfiilly submit that by the term ''axial dispersion 
coefficient," they do not seek any special or broader definition than that accorded to 
the term by one of skill in the art and that this is a term of art that is well substantiated 
as explained below. Use of this term in the Applicants' specification and claims is 
consistent with the ordinary and accepted meaning of this term. 

In order to determine whether or not a given process falls within the metes and 
bounds of the claims (specifically, whether a given process is conducted in a manner 
that yields a plug flow index of more than 50), one would only need to measure or 
model the axial dispersion coefficient for a particular implementation. A variety of 
techniques are known for measuring and modeling axial dispersion coefficient. As a 
illustration of the methods known in the art for determining axial dispersion 
coefficient applicable to the application, please see the following references and 
references cited therein (i) pages 529-531 and 532-551, Chemical Reactor Analysis 
and Design, 2"^ Edition, Gilbert Froment and Kenneth Bischoff, John Wiley & Sons, 
NY; ISBN 0-471-51044-0; (publication year: 1990); (ii) pages 2248-2250, Axial 
Dispersion in a Kureha Crystal Purifier, K Otawara and T Matsuoka, Journal of 
Crystal Growth, Volume 237-239 (publicafion year: 2002); (iii) pages 41-51, Reactor 
Models, S Dutta and R Gualy, Chemical Engineering Progress, (publication year: 
October 2000). While these references teach specific techniques, any method now 
known or later discovered may be used to quantify the axial dispersion coefficient 
which, in turn, can be used to determine whether a particular process is falls within 
the scope of the instant claims. 
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C. Conclusion 

In view of all of the above, claims 1-4, 6-9 and 20-26 are believed to be 
allowable and the case in condition for allowance which action is respectfully 
requested. The references that were cited and not relied upon are believed to be no 
more pertinent that those references that were relied upon. 

No fee is believed to be required by this response as determined on the 
accompanying transmittal letter. Should any other fee be required, please charge 
Deposit 50-1 123. Should any extension of time be required please consider this a 
petition therefore and charge the required fee to Deposit Account 50-1 123. 

Respectfully submitted, 



Date: November 3, 2003 




Stuart T. Langley #33,94Cky 
HOGAN & HARTSON LLP 
One Tabor Center 
1200 17th Street, Suite 1500 
Denver, Colorado 80202 
Phone: (720) 406-5335 
Fax: (720) 406-5301 
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product exits as crystal from the top of the 
column. Typical usage of a KCP is depicted in 
Fig. 3. 



3. Theoretical 

For scaling up the KCP» the following two 
phase counter current dispersion model [5] has 
been adopted as the steady state model. 



(1) 



The boundary condition is expressed as 

- ^ = VAC,,-. I. •- C,-/ ) at z = L, 



-r^ = 0 at 2 = 0, 



-Ey = Vy(CyB - C„....o) at r = 0, 



dCv 
az 



(2) 

(3) 

(4) 
(5) 

(6) 



where C, V\EAki^, a% and y are concentration of 
impurities, velocity, axial disp ersion coefficient, 
sweating rate constant, liquid, and solid, respcc- 
Lively. The model simplified various aspects of 
phenomena in the KCP: For example, temperature 
profile along the axis, recrystallization, si^e dis- 
tribution of crystals or washing effect is not taken 
into account. Currently, wc arc trying to develop 
models which incorporate more details of such 
phenomena than the present model. This will be 
reported in the future. 

The model has three parameters, i.e., sweating 
rate constant and axial dispersion coefficients for 
liquid and solid. Among these, axial dispersion 
coefficients can be estimated by well-known tracer 
techniqu e as discussed in the succeeding section. 
ihc sweating race constani has been estimated by^ 
analyzing impurity concentration profile along the 
axis. This will be published elsewhere. 



4. Experimental procedure 



To estimate axial dispersion coefficients and 
mean velocities of Hquid, tracer technique expe ri 
ments have been carried out for obtaining 
residence time distributions of liquid with the 
variatrons of tne size ot the KCP. 

When an idealized tracer pulse is introduced 
into a vessel with the boundary condition of the 
closed system, an axial dispersion coefficient can 
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Fig. 3. Representative flow tiiiigram of a KCP syslem. 
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be calculated by (see, e.g., [6]) 



(7) 



w here (Tq is the normalized sUtnclard deviation for 
the trac er response at t he exit of the vessel. 

Klonoclilorobenzene has been chosen as a 
tracer. Its temperature has been kept around the 
melting point of para-dichiorobenzene (/7-DCB) to 
minimize the disturbance lo the operation of the 
KCP. A certain amountof the tracer, e.g., 1 kg for 
the KCP of 0.51 m in diameter, has been fed as a 
pulse input from the top of the column after the 
KCP was stabilized to produce pure />-DCB. Then, 
the bottom product has been sampled every 3 min 
and its composition has been analyzed by a gas 
chromatography. To estimate the effect of the 
column diamclcr on the axial dispersion coeffi- 
cient, the similar experiments have been carried 
out for KCPs of 0,08 and 0.15 m in diameter as 
well with less amount of the tracer. 
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Hg, 4. Residence lime dislrihution of tracer as Uie liquid phase 
in the ICCP. 



5. Results and discussion 

Fig. 4 illuslrates the composition profile of the 
tracer of the sample taken from the bottom of a 
KCP whose diameter is 0.51 m. It appears that the 
flow of liquid phase in the KCP is not totally plug 
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Fig. 5. Comparison of axial disipcrsion cwirficients of liquid 
with « varitUion of inner di ameter of crystalJization columns. 

flow, but has a small degree of ihe back-mixing. 
The axial dispersion coefficient, i.e., £, of KCP 
with the variation of the column diameter, i.e., D, 
of the KCP is shown in Fig. 5. From this, the 
following relation between the size of the diameter 
of the KCP and the axial dispersion of the liquid 
phase in the KCP has been derived. 

E = kED'^. (8) 

In this expression, /ch is a proportional constant. 
This indicates that the larger the KCP is, the more 
significant the nxial dispersion is. Fig. 5 also shows 
thai ihc axial dispersion in the KCP can be 
substantially smaller than those reported for other 
types of column crystallizer for a given column 
diameter [7-8]. Thus, KCP might possibly oulper- 
form them . especially when their diameters are 
large. 

Experiments to investigate the axial dispersion 
coefficients of the solid phase have been executed 
for KCPs with the column diameter of 0.08 and 
0.15 m. Nevertheless, so far there has been no 
opportunity for such experiments using KCPs of 
industrial sizes, since soHd tracer, which arc 
proved to be not harmful to p-DCB and KCP 
processes, have not been found at this moment. 
The axial dispersion coefficients of the solid phase 
will he explored and reported in the future. 
Nevertheless, wc have carried out experiments 



with KCPs of 0.08 and 0.15m in diameter and 
have obtained rough effect of the column diameter 
on the axial dispersion coefficients in the solid 
phase. 

In designing the KCPs, it is most important to 
determine their column heights. With the model 
discussed in the theoretical session, if we estimate 
parameters included in the model by the experi- 
mental correlation equation obtained in Fig, 5 and 
the allowable concentration of impurities of the 
product is specified, the required heights can be 
calculated. The largest KCPs, specifically, 0.76m 
in column diameter, were constructed in 1998 in 
this manner. These KCPs have been successfully 
operated and their capacities for purifying ;?-DCB 
have been proved to be 13,000 T/y. Small axial 
dispersion coefficients have rendered possible to 
develop such KCPs. 



6. Concluding remarks 

The axial back-mixing in the equipment affects 
the performance of a column crystallizer. In this 
paper, the axial dispersion of the liquid phase in 
the KCP has been investigated and il is found how 
the diameter of the KCP affects the axial disper- i 
sion coefficient of the KCP. These coefficients have \ 
been compured to those obtained for other types j 
of column crystallizers and it has been found that i 
the former is substantially smaller than the latter. ■ 
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grams usually cover deflagration^ 
which is caused by sudden gas expan- 
sion by an accelerated reaction. The 
maximum pressure buildup due to such 
an explosion can be estimated from the 
pressure/'volume relationship of the 
complete reaction/'combustion process. 
The detonation process, on the otlier 
hand, is more violent, generating shock 
waves that travel at a speed several or- 
dej:s of magnitude higher than tlie pres- 
sure wave caused by deflagration. The 
local pressure buildup from this phe- 
nomenon, thus, is substantially higher 
— for' most hydrocaxbon/air mixtures, 
about 15-20 times the initial pressure, 
compared to a maximum of two Limes 
in case of deflagration. Sec, for exam- 
ple, Ref 6 for details on these explo- 
sion phcnomcnu. 

Therefore, the detonation charac- 
teristics of the reaction mixture 
should be studied, along with defla- 
gration, to address reactor safet}' 
more thoroughly. 

STEP 6 ~ DEFINE REAC- 
TOR TYPE AMD OTS HY- 
DRODYOyAMDCS. The reactor 
type is defined by: (a) the physical 
configuration of tlie volume occupied 
by the reaction system; (b) the flow 
mode of various streams in and out of 
the reactor; and (c) the hydrodynamic 
representation of the flows within tlie 
reactor volume. Thus, for example, a 
CSTR usually represents a stirred 
vessel with continuous flow of a ho- 
mogeneous fluid stream or continu- 
ous phase in and out of the vessel. 
The hydrodynamic representation of 
a CSTR a,ssumes the fluid stream is 
completely mixed as soon as it enters 
the vessel nnd attains the outlel com- 
position instantaneously. A PFR em- 
bodies the other extreme. It usually is 
tubular (with a high ratio of length, L, 
to diameter, D) in configuration. The 
flow in and out this reactor type is the 
same as that in the CSTR. The repre- 
sentation of the PFR hydrodynamics, 
however, assumes a total absence of 
axial mixing of the flowing streams 
within the reactor 

A homogeneous reactor usually is 



designed using the performance equa- 
tions of either the CSTR or PFR. This 
allows easier analysis, design, and 
scale-up of such reactors. But, these 
idealized reactors do not closely 
enough mirror the behavior of a real 
reactor Such a reactor is better repre- 
sented by a system with axial disper- 
sion or by an axially dispersed reactor 
(ADR), the performance of which 
falls between that of a CSTR and a 
PFR. For a simple first-order reac- 
tion, A B, the steady-state mass- 
balance equation of an ADR can be 
described as: 

Z>„^._4k_fe^ = 0 (1) 

where is the axial dispersion 
coefticienL 

The predicted behavior of the 
ADR approaches that of a PFR for no 
axial dispersion, /. e., = 0, and that 
of a CSTR for infinite dispersion 
(complete mixing), /.e., = oc. For a 
real reactor, D^^ is a finite number 
greater than zero. 

In developing a robust model, it 
may be advisable to start with an 
ADR equation, rather than a CSTR or 
PFR one. For initial model develop- 
ment, the ADR equation can be 
turned into simple CSTR or PFR ones 
by substituting the two limiting val- 
ues of D„. But, including the term 
allows you to fine-tune the model at a 
later stage to predict the reactor per- 
formance more precisely. To use this 
approach, however, the va hjc^nf /\ 
neeqs to be determincdT^v expgnT" 
ment or from one of the many corre- 
lations available in literature, for ex- 
ample, in Refs. 1^ a n d 7-IO> T hese 
correlations usually are available in 
tlie form of the dimensionless Peclet 
number, A';:^, (= Lu/D^,, where L and u 
represent tube length and fluid veloci- 
ty, respectively). The Peclet number 
typically is expressed as a function of 
Reynolds number, iV^^, which repre- 
sents the flow behavior within the re- 
actor. Nf^ > 100 normally indicates an 
approach to the behavior of a PFR, 
while N,^ < 1 to that of a CSTR, 
Qualitiitively, a longer tube length or 



higher fluid velocit>' means more 
PFR-Jike performance, and the re- 
verse more CSTR-like. 

The above discussions are for con- 
tinuous flow reactors. For both homo- 
geneous and heterogeneous reactors, 
however, part of the reaction mixture 
or one or more of the phases can be 
captive or in nonflow condition, 
while the rest can he in continuous or 
intermittent flow mode. Rcactants can 
be injected along the reactor length or 
reactor path in the case of multiple- 
reactor systems. The reactor product 
can be recycled with the reactor feed 
before or afl:er separation stages. 
These represent batch, semi-batch, 
multi-injection, and recycle reactor 
types. For heterogeneous reactors, the 
flow of the phases can be co- or 
coiuitcr-currcnt. 

For gas/solid reactors, the gas usu- 
ally is in a continuous flow mode. 
The solid, however, can remain fixed 
in position within the reactor through 
which the gas percolates. Or, botli 
solid and gas can flow continuously 
in CO- or counter-current mode 
through the reactor. With die solid re- 
maining fixed within the reactor, 
there arc again various reactor types 
— for example, the nonadiabatic 
packed-bed tubular reactor (NAPB- 
TR), adiabatic packed-bed multistage 
reactor (APBMSR), monolithic reac- 
tor, and radial flow reactor. An 
APBMSR can be operated widi direct 
quench (interstage cooling or heating) 
by feed or recycle-gas injection or in- 
direct quench by a heat exchanger. 
For continuous flow of botli gas and 
solid, the possible reactor types in- 
clude moving bed (co- and counter- 
current), bubbling and turbulent fiu- 
idized beds, circulating and fast flu- 
idized beds, and entrained beds. 

For gas/liquid systems, the reactor 
generally used is either a stirred tank 
or a bubble column. The bubble col- 
umn can be vertically or horizontally 
sparged, vertical or horizontal flow, 
gas-Iifled (internally or exlcmally), 
and with or without a forced liquid- 
circulation loop. The bubble columns 
also can be categorized according to 
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CataJytic cracking of pasoil {Continued) 
industrial performance of, 567, 592, 588 
kinetic model for, 20, 586 
kineUcs, 586 
in riser, 592 
Catalyhc cracking of //^hexiidecanc, (A 
Catalytic reaction: 
kinetic equalions for, 70 
mechanijims for, hi 
Catalytic roiictors: 
fixed bed, .194 
fluidized i>od, 567 
for kinetic analysis. 84 
Chain rcactit>ns, 26 
C'hain transfer, 33 
CTii-square test, 103 
Coke combustion, 209 
Coke deposition in chemical reactors: 
examples of, 442, 491 
modeling of. 437 
Coking, kinetics of: 
in bulcnc hydrogonation, 247, 25(1 
in catalytic cracking of gas oil, 243 
inside catalyst particle, 237, 239, 243 
in //-pentai^e isomerization, 244 
Complete segreg3tioi\, 526 
Complex reaction, activation, energy of, 3fi. Sec 
itho C(mseculive reactions; I'nrailcl 
reactions 
in gas liquid systems, 271 
kinetic equations for, 14 
Complex reaction networks, U» 
HJD method for, 45 
in perfectly mixed reactors. 360 
Woi-Pratcr treatment of, 16, 43 
Coniidencc region, 97 
Configiirational diffusitm, 1.S4, 
reaction accompanied by, 181 
Cxmsecutive reactions, W 
in batch reactors, 327 
with catalyst deactivation, 227, 443 446 
gas~lic]uid, 271 

in perfectly mixed reactors, 363 

with pore diffusion, 175 

in stirred tanks, 365 
Continuity equations: 

general' form of, 295 

simpJi/icd forms of, 298 
Continuous flow stirred tank renctor, 358 
Conversion, 5 

■and segregation, 527 
Cracking, sec Catalytic cracking;; I herfnal 
cracking 

Damkiihler number, 450 

for poisoning, 223 
Danckwerts' age distribution function, 275, 
l^anckwerls' boundary condition, 535 
Deactivaiicm oi catalysts, by coking, 230 

influence on autothermic operation, 428 

influence on selectivity, 243 

non-steady state roaclor behavior due to 
437, 491 

by poisoning, 220 
Deactivation constant, 232 
Deactivation function: 

for coking, 232 

for poisoning, 221 



Oecomposition for ester, batch reactor desicn 

lor, 312 ° 
Hegrce of advancement, 5 
Degree of polymerization, 33 
Dehydrogenation: 
of 1-butenc. sec Butrne dehvdrogen«tion 
of clhanol, stY Hlhjnol dehvdroi;eniifion 
Dcsjgn, sec Keactor dcsi^-n ' 
Differential method of kinetic analvsis, with 
batch data, 8, 42, 307 
for caiaiytic reactions. 87 
with tubular reactor diita, 339 
Differential reactor, 86 

Diffusional falsificaHon of kinetic parameters, 

10i5, Iu7 
Diffusion resistance: 
with bifunctional catalysis, 167 
in carbon to carbon -dioxidi- reaction 170 
with complex reactions. 175 
in gas-liquid reactions, 25ft 
in gas -solid reactions: 
**xtcrnal, 132. 171 
internal, see Pore diffusion 
Ditftisivity, siv Liffective diffusivity 
Discrimination, among rival models. KH 
Dispersion coefficient, 530 

determination of. 530, 533, 547 
Diverpna! rrilerion. in sequential planninc 
104 ^' 
Drag coefficienl, in fluidi/ed beds, 570 
Driving force groups in. catalytic reactions, 79 
l->ual function catalysts, >ee Bifuncliniial 
catalysts 

Dynamic programming, applicalion to 
mullibed adiabatic reactor design. 419 

liffeclive binary diffusivity, 130 
f^ffeclivi* diffusion model: 
in fixt?d bed reactors, 451 
in multiphase flow reactors, 608 
nffective diffusiviiy: 
in catalyst particle, 145 
in packed beds: 
in axial direction, 447 
in radial direction, 45 1 
in turbulent jlow, 297 
fcifioctivcness factor, 157, 478 

global, 171, 186 
LYfeclive surface diffusiviiy. 154 
nffective thermal conductivilv: 
inside catalyst pellets, 183' 
in packed beds: 

in axial direction, 448 
in radial direction, 452 
fc'nerg)' of activation, see 

Activation energy 
Knergy equation: 
geneial form of. .301 
simplified forms of, 302 

Rnhancement factor for gas absoq^tion, 262, 
266 

tiolvos number. 622 

based on bubble diameter, 637 
tquivalent i-eactor volume, 341 
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activation energy, 36 
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MOOELS ACCOUNTING FOR TLOW PAITFRNS 

Similar approaches have been taken by Taiong and Methot (1976). Also, 
Danckwerts (1958) defined an index of segregation based on the variances of 
the ages in the fluid elements and the vessel as a whole; however, there is nu 
direct relation with reactor performance. 

Again, no general correlations of these parameters are available. Surveys 
of various aspects of relating turbulent mixing concepts to micromixing in re- 
actors appear in Brodkey (1975). Surveys of the topic of mixing and chemical 
reactions are given by Olson and Stout (1966) and by Villermaux (1976a, 1986b). 

A final important point is that all of the above is valid only for isothermal 
conditions in single-phase systems. Extension to other cases requires the intro- 
duction of interactions with the second phase and/or heat exchange walls, for 
example, and the age distribution and micromixing functions depend much more 
on the details of the system. A formal treatment would use joint probability 
distribution functions, but this rapidly gets extremely complex. The population 
balance models can give some insight into the two-phase situation. They are 
discussed below. 



1 2«6 

MODELS ACCOUNTING 
FOR FLOW PATTERNS 
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12.6.1 Axial Dispersion Model 

and Tanks-in-Series Model 

The above methods are quite general in that they are not based on any specific 
physical model. On the other hand, they are limited in presuming that the RTD 
is available and also by restriction to homogeneous phases. Thus, other tech- 
niques are required for predictive purposes and design where the actual reactor 
is obviously not available for a tracer test. To develop general correlations of 
behavior, the usual engineering tactics of utilizing mathematical models with 
parameters to be determined from the experimental data are useful. These pa- 
rameters are then correlated as functions of fluid and flow properties, reactor 
configurations, along with other important features, and can be used in design 
calculations. The flow models are senuempirical in nature, but hopefully will 
bear some relation to the actual flow patterns in the vessel. 

There are several types of models that have found useful applications, but 
two are most common — at least when found to be adequate to represent the 
physical situation. The first is usually termed the "axial dispersion" or "axial 
dispersed plug flow" model [Lcvenspiel and Bischoff (1963)1 and takes the form 
of the one-dimensional diffusion equation with a convective term. The second 
is a series of perfectly mixed vessels, some features of which have been discussed 
in Sec. 12.3. 

The mass balance for the axial dispersion model is, from Chapter 7, 



+ w-r = :- ~ + r{c) 
dt flz c*z f)z 



(12.6.1-1) 



In Eq. 12.6.1-1, u is taken to be the mean (plug flow) velocity through the vessel, 
and is a mixing-dispersion coefficient to be found from experiments with the 
system of interest. One important application is to fixed beds, as discussed in 
detail in Chapter 11, and then it is usually termed an "effective" transport model, 
with Da - Dca- However, the axial dispersion model can also be used to ap- 
proximately describe a variety of other reactors. 
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One of the main benefits of this model is its analogy to the diffusion 
equation, and the possibility of utilizing all of the classical mathematical solutions 
that are available [e.g., Carslaw and Jaeger (1959) and Crank (1956)J. 

One common approach to determinir^g t^ ^ modp) parameter [n 
perform a residence time distribution test on the reactor and to choose the 
value of Dfl so that the model solution and experimental output curve "afiree^ 
[ e>g.,^by least squares tec hniq ues (see Sec. 12.6.3)]. Figure 12.6.1-1 shows 
the t(e')-6' curves o^ the mode! tor an impulse input with "closed" boundaries 
(here, 6' - QF'/V = Gm/L, with L = length of the reactor): 



with a, positive roots of 



4a? 4re« + Pe^ 



tan a, = 



4Pe.a, 
4a, - PcS 



or 



where 



Pe., - uUD,. 




Fi9Hr« 12.6.1.1 

E curves in closed vessels (or various exlenls of dispersion; Pe, = uUD». Attor 
Levenspi'el (1 972). 
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he diffusion 
cal solutions 
)]. 

ter Da is to 
choose the 
curve agree 
6.1-1 shows 
boundfiries 



(12.6.1-2) 



C 
M/V 



(See Carslaw and jaeger (1956).] A thorough evaluation was given by Brenner 
(1962). It is seen that the axial dispersion model can represent mixing behavior 
ranging from perfect mixing (U, ^) to plug flow (D,, 0). 

Other t)'pes of boundary conditions can also be used when solving Hq. 
12.6.1-1 (some details are given below when discussing chemical reactor appli- 
cations); for example, in an "infinite" pipe, the dispersion characteristics up- 
stream from the tracer injection would be identical to those within the system, 
and a different mathematical solution results: 

The £ curves from Eq. 12.6.1-3 arc similar to those from Eq. 12.6.1-2 for rea- 
sonably large values of Pe«. 

Many of the available solutions are given in Levcnspiel and Bischofi (1963), 
Himmelblau and Bischoff (1968), and Wen and Fan (1973). As a practical matter, 
these mathematical complexities should probably not be taken too iilerally be- 
cause the precise conditions at the boundaries of real equipment cannot usually 
be exactly defined, in any event. That is, if the boundary condition details give 
significantly different results, be cautious about the axial dispersion model as 
being a faithful representation of the flow patterns in the vessel. These differ- 
ences are only seen for large values of D,, (actually, small PeJ, but the physical 
basis is a diffusionlike situation, which would not be valid for extremely large 
levels of mixing. 

EXAMPLE 12.6.1.1 

AXIAL DISPERSION f>/ODEL hOR LAMiNAR FL OW IN ROUND TUBES 

This is one situation where the exact convective diffusion equation can be solved 
and compared with the one-dimensional dispersion equation. The complete 
equation is 



(a) 



where u is the mean velocity. 

Taylor (1953) has discussed this in detail, and a summary of the results is 
given here. Taylor and Aris (1956) showed that the main variable of interest, 
the mean concentration, is 



2 



rC dr 



could be found by averaging Eq, (a) (for impermeable walls): 



a(C) 
a/ 



#(C> 



with 
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(b) 

«:) 

(d) 



See Hunt (1977) for a more recent exposition. The step from Eq. (b) to Hq. (c) 
is only possible under conditions that the radial concentration proHle has a 
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